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ABSTRACT 

We explore the detectability of bulk motions in the X-ray emitting intracluster medium (ICM) using 
a catalog of 1,836 mock Astro-E2 observations of simulated clusters of galaxies. We generate high 
resolution mock spectra for two observing strategies: a four-pointing mosaic and a single central pointing. 
Normalizing to 200 (400) photons in the iron Ka region for the mosaic (central) study, we fit Poisson 
realizations of each simulated spectrum to a velocity broadened isothermal plasma emission model. We 
confirm that the velocity characteristics (mean and dispersion) returned by the spectral fittings are 
unbiased measures of the emission-weighted values within the observed region, with scatter ±55 km s“^. 
The maximum velocity difference between mosaic element pairs Aumax has ^ 6% likelihood of being 
transonic (Aumax > 0.5cs), and the likelihood falls steeply, p oc (Aumax/cs)~"^, at high Mach number. 
The velocity broadening parameter ay from the central pointing fit exceeds the thermal value in 49% of 


the cases, with again a cr„ tail at large dispersion, 
the strongest signal for each observing strategy. 

Subject headings: clusters: general — clusters: ICM 
medium — X-rays: clusters 

1. INTRODUCTION 

In a hierarchical model of structure formation, galaxy 
clusters are built from the continuous accretion of sub¬ 
systems and, therefore, clusters are expected to frequently 
display visible signatures of mergers. Binary clusters 
(Jones & Forman 1984) offer evidence of precursor sys¬ 
tems while clusters with multiple emissivity peaks, strong 
isophotal twists and ellipticity variations possess mor¬ 
phologies associated with the later stages of merger evo¬ 
lution (Mohr, Fabricant & Geller 1993; Buote & Tsai 
1996). Bulk motions and turbulence of the hot intraclus¬ 
ter medium (ICM) are thus expected to be evident in line 
emission from the plasma, in particular from helium-like 
iron, FeXXV (Dupke & Bregman 2001a,b; Sunyaev et al. 
2003; Inogamov & Sunyaev 2003). 

X-ray satellite missions to date have had spectrometers 
with energy resolutions marginally adequate for the ob¬ 
servation of cluster velocity structure (Dupke & Bregman 
2001a,b). The new microcalorimeter technology employed 
in the X-ray Spectrometer (XRS) instrument aboard the 
Astro-E2 satellite, however, is designed to have an energy 
resolution of 6.5 eV FWHM across its entire frequency 
range (Furusho et al. 2004), giving a response contribu¬ 
tion to the line broadening of ~1.5 eV and finding line 
centroids to within about 2 eV. This implies that the in¬ 
strument will be sensitive to velocity structure down to 
~100 km/s near the Fe Ka line. Since the infall velocity 
of material accreting onto a large cluster can be an order 
of magnitude larger than this limit, Astro-E2 will enable 
unambiguous observation of gas motion in dynamically ac¬ 
tive systems. 

Cosmological simulations provide a framework for inves¬ 
tigating the statistical properties of the internal gas veloc- 
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ity distribution. In addition, access to the complete 3-D 
structure and time evolution of the simulated clusters can, 
through comparative studies of observations and simula¬ 
tions, help unfold the evolutionary history and dynamical 
state of observed clusters. 

Until now, no estimate of the frequency with which clus¬ 
ters should show detectable velocity structure has been 
made, although Fujita et al. (2005) have used mock obser¬ 
vations of simulated clusters to evaluate the detectability 
of turbulence in a specific case study. In this paper, we 
evaluate the fraction of clusters expected to show observ¬ 
able structure in the Fe Ka line by performing mock Astro- 
E2 observations of a catalog of simulated clusters. The 
simulations model clusters in a ACDM cosmology with a 
preheated assumption for the ICM that matches observed 
local scaling laws (Bialek, Evrard & Mohr 2001). 

In §2, we describe the simulation sample and the strate¬ 
gies we use to produce mock observations. Results in the 
form of statistical expectations and illustrative examples 
are presented in §3. We summarize our findings in §4. 

2. MOCK OBSERVATIONS 
2.1. Preheated Cluster Simulations 

We use a catalog of 68 clusters evolved with P3MSPH 
(Evrard 1988) under a preheated assumption for the ICM. 
The level of entropy introduced into the initial conditions 
is tuned to 105.9 keV cm^ in order to match the observed 
scalings of luminosity and ICM mass with temperature 
(Bialek et al. 2001). 

The simulations were produced using a multi-step pro¬ 
cedure outlined in Bialek et al. (2001). The underlying 
cosmology is a flat, concordance model with flm = 0.3, 
Da = 0.7, Df, = 0.03, as = 1.0, and h = 0.7, where the 
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Hubble constant is defined as 100ft. km s“^ Mpc“^ and ag, 
is the power spectrum normalization on 8ft“^ Mpc scales. 
Details of the full ensemble will be presented in Bialek, 
Evrard & Mohr (2005). One member of the ensemble con¬ 
tained a prominent cold front feature, generated during a 
merger by the separation and resultant adiabatic cooling 
of the core of an infalling satellite, as discussed by Bialek, 
Evrard & Mohr (2002). 

The cosmological values assumed for the models are in 
line with WMAP allowed parameters, with the exception 
of the baryon fraction. The value flb/rtm = 0.1 used in 
the simulations is lower than the WMAP value 0.17 ± 0.01 
(Bennett et al. 2003). The mass fraction in intracluster 
gas is expected to be less than the cosmic baryon frac¬ 
tion due to galaxy formation and energy exchange with 
dark matter during mergers (Thomas & Couchman 1992). 
The former process removes ~ 15 — 20% of the baryons 
from the hot phase in rich clusters while the latter expels 
~10% of the gas from the potential (Erenk et al. 1999). 
Still, the model gas fractions are likely to be somewhat low 
compared to current observational estimates. However, in 
the results presented here, the baryon fraction affects only 
the normalization of the ICM X-ray spectrum. Our mock 
spectra are tuned to a fixed line normalization, which at 
fixed cluster mass scales oc (Hb/Dm)^texp, with texp the ex¬ 
posure time. A higher or lower baryon fraction can thus 
be absorbed by appropriate rescaling of the exposure. 

The configuration of each simulation is stored at twenty 
output times, spaced equally in time from the initial red- 
shift Zi = 20.82 to the present. We employ the final nine 
outputs of the 68 models in this study, corresponding to 
the redshifts 0.540, 0.448, 0.365, 0.290, 0.222, 0.160, 0.102, 
0.049 and 0.0. This yields 612 cluster realizations which 
we treat, under an ergodic hypothesis, as statistically in¬ 
dependent. The time separation of ~ 0.7 Gyr between 
outputs typically exceeds the crossing time at r 2 oo, justi¬ 
fying the ergodic assumption. 

We further enhance our sample by considering the line 
of sight (LOS) velocity structure of each realization along 
three perpendicular axes. Although these projections are 
linked by dynamics, the LOS velocity structure is gen¬ 
uinely independent. We thus have three projections of 
each of the 612 realizations, for a total of 1,836 exposures. 
The members of the simulated sample range in spectral 
temperature from 1.5 keV to about 8 keV, with cluster 
masses M 200 ranging from 0.015-2.4 xIO^^Mq. 

The derived spectra and images, along with associated 
parameters describing the simulations are publicly avail¬ 
able as part of VCE, the Virtual Cluster Exploratory^. 

2.2. Creating the Mock Exposures 

We create mock Astro-E2 spectra of these realizations 
following two observing programs illustrated in Figure 1. 
The first program involves a search for velocity gradients 
by measuring differences in the mean velocity of the Fe 
Ka lines for multiple pointings of the Astro-E2 instru¬ 
ment. This strategy allows one to measure both the magni¬ 
tude and a crude direction of the velocity gradient. Mul¬ 
tiple pointings are required because the average angular 
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resolution of the X-ray telescopes on Astro-E2 is ~1.9' 
(half-power diameter®) and the detector configuration is 
roughly a square with 2.9' on the side. Since less than 60% 
of the encircled energy is contained within 2', spatially re¬ 
solved spectroscopy using sub-regions of the detector will 
be highly contaminated with photons from neighboring re¬ 
gions. Therefore, to carry out this observational strategy, 
we create four mock spectra of each cluster projection in 
a 2x2 box centered on the cluster. The second program 
involves searching for extra-thermal broadening of the Fe 
Ka line using a single, central pointing. 



map for a typical cluster in our catalog. The black boxes indicate 
the fields for the four-pointing mosaic program while the white box 
shows the central field of the line broadening program. Each field 
is 2.9 arcmin on a side, the FOV of Astro-E2, which corresponds to 
321 kpc at 2 : = 0.1. 

The thermal plasma emission spectra for this study are 
compiled using the APEC code (Smith et al. 2001) from 
the XSPEC 11.3.1 suite of spectral analysis tools (Arnaud 
1996). APEC allows us to include thermal broadening of 
the line spectra by setting the APECTHERMAL toggle 
to ‘yes’. Further, XSPEC allows us to fold in the antici¬ 
pated XRS response functions®. We use XSPEC/APEC 
to write a reference table of spectra in flux units at 176 
temperatures spaced logarithmically to cover our range of 
interest. This spectral table is limited to an energy range 
appropriate for studying the Fe K complex. The final spec¬ 
tra are limited to the interval 5.9 keV to 6.4 keV. All spec¬ 
tra are generated with the cluster at a fiducial redshift of 
0.1, so that the FeXXV Ka line is centered at 6.09 keV 
and the FeXXVI Ka line at 6.33 keV (the hydrogen-like 
iron line is actually two separate lines with centers at 6.32 
and 6.34 keV as a result of fine structure (Verner, Verner 
& Ferland 1996)). Figure 2 shows the Ka line region of 
reference spectra at a few relevant temperatures. 

To create a simulated spectrum for a mock observa¬ 
tion, we interpolate on the reference flux table to generate 
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an emission spectrum for each gas particle of the simula¬ 
tion contained in the Astro-E2 FOV. There are typically 
~ 1000 particles in the the FOV, and their contributions 
are summed to obtain the complete spectrum. To convert 
the resulting flux spectrum into an XRS count spectrum, 
we first define an exposure time that requires either 200 or 
400 total counts in the Fe Ka line region within the energy 
range of interest, i.e., 6.0-6.2 keV. (The 200 count crite¬ 
rion applies to the four-pointing mosaic study, while 400 
counts are used for the central pointing spectra.) A sin¬ 
gle, discrete Poisson realization of each spectrum is then 
created, using the poidev routine of Press et al. (1992). 
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Fig. 2.— APEC thermal spectra of 0.4 solar metallicity plasma 
(arbitrary flux units) in a narrow energy range centered on the iron 
Kq: complex at redshift 2 = 0.1 are shown for typical temperatures 
of plasma in the simulations. The Astro-E2 XRS energy response 
has been folded in. 

2.3. Additional Assumptions 

The iron abundance of the ICM is not determined by 
cluster simulations, and so has to be assumed based on 
data from observed clusters. Given the small field of view 
of the XRS on-board Astro-E2, we are probing regions 
fairly close to the cluster’s core (< 0.2 r 2 oo)- The iron 
abundance in cluster cores is typically higher than the av¬ 
erage over the whole cluster, due to the frequent presence 
of central metal abundance gradients in cold core (“cool¬ 
ing flow”) clusters (Ulmer et al. 1987; White et al. 1994), 
which constitute the majority of clusters. Assuming an 
average central iron abundance of 0.27 solar for non-cold 
core clusters and 0.47 for cold core systems (DeGrandi et 
al. 2004), and recognizing that 70 — 90% of an X-ray flux- 
limited sample of clusters have cold cores (Edge, Stewart & 
Fabian 1992), we obtain an average central iron abundance 
of ~0.4 relative to the solar photospheric values of Anders 
& Grevesse (1989). This is the value assumed throughout 
this work. 

One of the main obstacles in measuring velocities of 
the intracluster gas with current spectrometers is the in¬ 
ability to accurately calibrate the temporal and spatial 
fluctuations of instrumental gain, the conversion of pulse- 
height into photon energy (Dupke & Bregman 2001b). 


For the case of CCDs the intrachip (positional dependent) 
gain changes are related to the charge transfer inefficiency, 
which evolves with time. Often these fluctuations are on 
the same order as the velocities one is trying to measure. 
Although this is not an issue for the XRS calorimeter, 
small drifts in the temperature of the detector heat sink 
can cause global gain variations. The absolute precision of 
the gain variations in the XRS is expected to be ~ 1 — 2 eV 
(50—100 km s“^) at 6 keV and will have to be inserted into 
the uncertainties of the velocity measurements described 
here (Figueroa, personal communication, 2004). These un¬ 
certainties were not incorporated in our evaluations of ve¬ 
locity gradients, but they do not affect the overall results 
derived in this work. Similarly, for the purpose of this 
paper we do not fake a background spectrum. The main 
source of background for the XRS is due to energy deposi¬ 
tion in the detector by energetic protons. However, given 
the small detector size and the distribution of the proton 
energy spectrum, its contribution in the frequencies of in¬ 
terest (< 10 keV) are expected to be negligible (< 2-5%) 
for our mock observations. 

Finally, in this paper we assume that the intracluster 
plasma is optically thin. Gil’fanov, Sunyaev & Churazov 
(1987) have pointed out that this may not always be true 
for the Fe Ka line. Resonant scattering (the absorption 
and immediate re-emission of an Fe Ka line photon by an 
Fe ion) can be significant in the cores of clusters. If so, the 
line emission from the cluster core is reduced because of 
the scattering of photons out of the field of view. Early ev¬ 
idence for resonant scattering (Molendi et al. 1998) seen in 
the Perseus cluster has proven ambiguous, allowing other 
possible interpretations such as an overabundance of Ni 
due the central SN la ejecta dominance (Dupke & Arnaud 
2001). XMM observations of the Perseus cluster have not 
exhibited resonant scattering in the Fe Ka line, and the 
absence of observed scattering has been used as evidence 
for the presence of turbulent (non-thermal) velocities in 
clusters’ cores (Gastaldello & Molendi 2004, Ghurazov et 
al. 2004). 



Fig. 3.— Histogram of the optical depth for the central pointing 
of the 1836 cluster projections used in this study. 
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To quantify the potential impact of resonant scattering, 
we compute the optical depth for the central pointing of 
each cluster projection used in this study (see Appendix 
A for details). To be conservative, we calculate the optical 
depth through the entire cluster. The results are shown 
in Figure 3. When the non-thermal velocities seen in the 
emission line are included, the average optical depth for 
the FeXXV Ka line in the FOV of the central pointing is 
0.59, with standard deviation 0.34. The results of Chura- 
zov et al. (2004) and Gil’fanov et al. (1987) show that 
for r ~ 0.6, the flux and the linewidth near the cluster 
core will show approximately a 15% reduction. Addition¬ 
ally, it is possible that some deformation of the line profile 
may occur (Gil’fanov et al. 1987). None of these effects 
would introduce errors in the line centroids we recover. 

The reduction in linewidth would affect our line broaden¬ 
ing study, but, as our results will show, a 15% reduction 
is not prohibitive. 

2.4. Velocity Gradients Via Multiple Pointings 

We first investigate the LOS velocity structure of the 
simulated clusters by creating one set of four mock expo¬ 
sures tiled around the most bound position of the relevant 
cluster, as shown by boxes 1-4 in Figure 1. The mock expo¬ 
sures result in four separate Poisson-noise-added spectra. 

We rebin these spectra to ensure at least 15 counts in each 
bin using the grppha routine from the FTOOLS package^. 

We fit the grouped spectra to find the position of the 
centroid velocity Vi of the Fe Ka line in image i. The 
spectra are fit in XSPEG to a BAPEG model. BAPEG 
is a thermally broadened APEG that allows for additional 
Gaussian velocity broadening under a parameter 

The resulting centroid positions are differenced to de¬ 
fine velocity gradients Vij = \vi — Vj \ among the six distinct 
pairs of the four pointings. Thus, from our 1,836 indepen¬ 
dent cluster projections we generate 7,344 spectra which 
are used to give 11,016 velocity differences. 

An advantage of working with simulated clusters is our 
complete knowledge of the underlying velocity structure 
of the gas. In an isothermal cluster, the redshift of the 
line center could be recovered from the density-squared 
weighted velocity (henceforth called the EM-weighted ve¬ 
locity) of the emitting gas. Our simulated clusters are not 
isothermal, but it happens that the Fe Ka line flux is not 
strongly temperature dependent within the temperature 
range expected for hot IGM gas. The peak flux varies by 
only a factor of 2 over a temperature range from 3 keV 
to 12 keV. For this reason, an EM-weighted velocity can 
be regarded as a predictor of the line center that would be 
recovered from a spectral analysis. Therefore, as a point of 
comparison we also generate the EM-weighted radial ve¬ 
locity of each Astro-E2 FOV directly from the simulation. 

To minimize systematic issues, we find the velocity of 
each FOV for a given cluster projection relative to the 
mean value of the four FOV’s for that projection. We re¬ 
peat this procedure using the BAPEG redshift values. A 
histogram of the differences between the velocities recov¬ 
ered in this way from the BAPEG spectral fits and from 
the EM-weighted averages is shown in Figure 4. Note that 
the XRS instrument is expected to find line centroids with 

^ see http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftools_menu.html 
® see ftp://legacy.gsfc.nasa.gov/astroe2/nra_info/astroe2_td.pdf 


a limiting resolution of about 2 eV near the Fe Ka com¬ 
plex®, which corresponds to an expected limiting accuracy 
of ~100 km/s near the Fe Ka line. The histogram of differ¬ 
ences has a dispersion of 54 km s“^. We conclude, there¬ 
fore, that 200 counts in the line (with grouping applied) 
results in velocity measurement accuracies very compara¬ 
ble to the expected 2 eV pre-launch instrument response. 
This result is consistent with the findings of Fujita et al. 
(2005). 



'^BAPEC '^EM (km/S) 

Fig. 4.— Histogram of the differences between velocities recov¬ 
ered from BAPEC fittings and those generated using EM-weighted 
velocities. 

2.5. Nonthermal Velocities Via Line Broadening 

Line broadening is another means of recovering the ve¬ 
locity structure of the emitting plasma. As described 
by Sunyaev et al. (2003) and by Inogamov & Sunyaev 
(2003), the large atomic mass of iron suppresses thermal 
line broadening, making the Fe K line width a sensitive 
probe of bulk gas motion. For line broadening studies, 
only one mock Astro-E2 image of each cluster centered on 
the most bound position is used, giving a total of 1,836 ob¬ 
servations. For this study, our final spectra are generated 
by requiring 400 counts in the line region (6.0-6.2 keV). 
The velocity broadenings of the spectra are returned as 
the parameter of the BAPEC fit. 

Once again, we perform a theoretical check on these re¬ 
sults. The EM-weighted velocity dispersion of the particles 
within the Astro-E2 field of view is taken to be a predictor 
of the best-fit BAPEC cr„ parameter. The differences be¬ 
tween these two measures of velocity broadening are shown 
in Figure 5. The dispersion for the entire data set of 1,836 
broadenings is 55 km s“^. 

The data of Figure 5 show evidence of a systematic 
effect. The mean value of ay{BAPEC) — av{EM) is 
clearly offset from zero. For the complete sample, it is 
—26.6 km s“^ (median —25.4 km s“^). This systematic 
is at least partly explained by the fitting process. The 
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spectral resolution (FWHM) of the XRS is about 6.5 eV, 
corresponding to a response contribution to the broaden¬ 
ing of Or = 1.5 eV or ~ 75 km s“^ near the Fe Ka line. 
This implies that XSPEC should be unable to recognize 
physical broadenings significantly less than 75 km s“^ in 
our spectra. Clusters with EM-weighted broadenings sig¬ 
nificantly less than 75 km/s are therefore reported as zero 
broadening in the BAPEC fit. 



odBAPEC)-odEM) (km/s) 

Fig. 5. — Histogram of the differences between the velocity broad¬ 
ening parameter returned by BAPEC and the EM-weighted velocity 
dispersion. The dashed line denotes the portion of the sample with 
BAPEC uncertainties > 5cr^, (BAPEC) (the “zero broadening” sam- 
pie). 

An investigation of this interpretation can be made us¬ 
ing the BAPEC-reported uncertainty in the broadening 
value returned. If we divide the uncertainty in the broad¬ 
ening parameter by the value of the parameter, we ex¬ 
pect to obtain a large number for spectra with broaden¬ 
ings that fall below the resolution of the XRS instrument. 
As a test, we construct this ratio and consider any spec¬ 
trum with a broadening uncertainty greater than five times 
the broadening value to belong to the “zero-broadening 
sample”. Implementing this definition yields a zero¬ 
broadening sample of 492 clusters which has a mean veloc¬ 
ity broadening reported by BAPEC of 9.0 x 10“^ km s“^ 
(median 1.8 x 10“"^ km s“^) and a standard deviation of 
0.025 km s“^. These numbers are well below the resolu¬ 
tion of the XRS instrument, confirming our intuition that 
a significant fraction of the overall cluster sample have 
physical broadenings unresolvable by XRS. The contri¬ 
bution of this zero-broadening sample to the histogram 
of Figure 5 is shown as a long-dashed line. Removing 
projections with zero reported broadening shifts the mean 
value of av{BAPEC) — ay{EM) to —13.4 km s“^ (median 
— 10.4 km s“^) and leaves the dispersion unchanged. 


3. RESULTS 

3.1. Velocity Gradients 


Figure 6 shows the cumulative frequency distribution of 
velocity differences, normalized to the sound speed in the 
ICM, that result from the four-pointing tiling program. 
Because the ICM temperature is determined by the grav¬ 
itational potential (j) of a cluster, the sound speed in the 
ICM gas is approximately 


Cs = 


15kTicM 




( 1 ) 


y Sfirrip 

where fiirip is the average mass per ICM particle. Since 
the bulk velocity structure of the ICM is determined by 
mergers also driven by the gravitational potential of the 
dominant cluster, one expects infall velocities 

Uinf « (2) 

All explicit dependence on (/>, hence on the size of the clus¬ 
ter, vanishes when infall velocity is scaled by the sound 
speed. 
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Fig. 6. — Cumulative likelihoods of the normalized velocity gra¬ 
dients Av/cs are shown for the case of only maximum velocity dif¬ 
ferences (solid) and the case of ah pointing pairs (dotted). A fit to 
the former is given by equation (3) (dashed). 

For each simulation, we determine the sound speed from 
the best-fit spectral temperature of an isothermal fit to the 
0.5 — 9.5 keV emission of plasma within rsoo, the radius 
within which the mean interior density is 500 times the 
critical value. This radius generally encloses several thou¬ 
sand of the gas particles making up the simulated cluster. 
The 0.5 — 9.5 keV spectra are created by interpolating on 
a reference table of isothermal plasma emission spectra to 
find the contribution from each gas particle, and then sum¬ 
ming these contributions to obtain the complete spectrum. 
The procedure is the same as that outlined in §2.2 for cre¬ 
ating mock Fe Ka line spectra, except that the table of 
reference spectra for interpolation is generated using the 
MEKAL plasma emission code (Mewe, Kaastra & Liedahl 
1995) and the spectral resolution is 150 eV per bin. The 
best-fit spectral temperature is determined by minimizing 
the chi-square of the residuals generated by differencing 
the complete 0.5 — 9.5 keV spectrum of the simulation 
with isothermal plasma spectra created using MEKAL. 

The solid line in Figure 6 shows the cumulative fre¬ 
quency distribution of the maximum velocity gradient 
among the six pairs of each four-pointing mosaic. The 
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fraction of cluster projections in the sample of 1,836 with a 
maximum normalized velocity difference larger than some 
value X is well fit by 

f{AVmax/Cs > x) « ^1 + (^) ^ • (3) 

The number of cluster projections in our sample with a 
maximum normalized velocity difference larger than half 
the sound speed for the cluster is 111, or approximately 
6.0% of the sample. The total number of splittings out 
of the 11,016 measured to be larger than O.Sc^ is 200, or 
approximately 1.8%. 

The likelihood of high Mach number splittings falls off 
dramatically, as {Av/Cs)~'^■ The lack of very high Mach 
collisions is expected from the arguments given in equa¬ 
tions (1) and (2) above. One anticipates Vinf ^ a few xc^ 
because larger Mach numbers could not be generated by 
hierarchical clustering driven by gravity. Our mock spec¬ 
tral measurements are consistent with this expectation. 
Note that the quoted value of the sound speed always refers 
to the dominant member of a merging pair. With respect 
to Cs of the smaller member, high Mach number collisions 
are certainly possible. But the measured “broad-beam” 
temperature, and sound speed, of the merging system will 
always be driven by the larger member, and thus the data 
presented in Figure 6 are appropriate for comparison with 
observation. 



0.01 0.10 1.00 
AV/C3 


Fig. 7.— Cumulative likelihoods of the normalized velocity gra¬ 
dients Au/cs are shown for the case of maximum velocity in three 
subsamples grouped according to simulation output redshift. The 
short-dashed line is the subsample of clusters at z > 0.365, the long- 
dashed line is clusters at 0.365 > 2 > 0.102, and the solid line is 
clusters at z < 0.102. 

It is important to recognize that even though all clus¬ 
ters have been moved to a fiducial redshift of z = 0.1 to 
construct mock observations, we have used clusters recov¬ 
ered from a simulated sample at several different output 
redshifts ranging from z = 0.540 to z = 0. It is therefore 
natural to ask if the frequency of large velocity gradients 
depends on the epoch of observation of the cluster sam¬ 
ple. To address this question, we have constructed the cu¬ 
mulative frequency distribution of the maximum observed 


normalized velocity gradient for three subsamples grouped 
by simulation output redshift. The results are shown in 
Figure 7. The Kolmogorov-Smirnov (K-S) test (see, e.ff., 
Press et al. 1992) returns a 48.9% likelihood that the high 
redshift and middle redshift subsamples follow the same 
distribution, a 44.5% likelihood that the low and middle 
redshift subsamples follow the same distribution, and an 
8.1% likelihood that the high and low redshift subsamples 
follow the same distribution. 

Another question to consider is whether the size of the 
Astro-E2 FOV relative to the cluster’s size plays a role 
in the frequency of large velocity splittings observed. We 
address this by constructing the cumulative frequency dis¬ 
tribution of the maximum observed normalized velocity 
gradient for two subsamples grouped by the physical size 
of the cluster (recall that all clusters are placed at the same 
fiducial redshift for observation, so that physical size cor¬ 
responds directly to angular size). Clusters with an r 2 oo 
larger than 1.46 Mpc (the median value for our sample) 
make up the large cluster subsample, while clusters with 
'C 200 < 1-46 Mpc form the small cluster subsample. The 
resulting distributions are shown in Figure 8. The K-S 
test gives a 58.6% likelihood that both subsamples follow 
the same distribution. 



0.01 0.10 1.00 
AV/C3 


Fig. 8.— Cumulative likelihoods of the normalized velocity gradi¬ 
ents Av/cs are shown for the case of maximum velocity in two equal 
subsamples grouped according to cluster size. The dashed line is the 
subsample of clusters with r200 < 1-46 Mpc, the solid line is clusters 
with r2oo > 1-46 Mpc. 

For purposes of illustration, we highlight the cluster 
which exhibits the maximum velocity gradient: cluster 
number 2 at output redshift 0.222 in y-axis projection. 
This projection yields a maximum normalized velocity gra¬ 
dient of 1.47 Cs- Data relevant to that cluster are compiled 
in Figures 9-13. Figure 9 shows Poisson realizations of the 
four-point mosaic spectra created using the procedure out¬ 
lined in §2.2 and grouped using the procedure explained 
in §2.4. These spectra show that the top two quadrants of 
our observing pattern are blueshifted relative to the bot¬ 
tom two. This accurately represents the underlying mo¬ 
tion of the cluster, as can be seen from Figure 10 which 
shows the perfect flux spectra generated from the Simula- 
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tion overlaid with the BAPEC fits to the spectra of Figure 

9. 
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Fig. 9.— The integrated and grouped Poisson-noised spectra used 
in the velocity gradient study for cluster 2 at 2 = 0.222, viewed along 
the y-axis. The spectra are generated with the cluster at a fiducial 
redshift of 0.1. Each panel displays the spectrum of the correspond¬ 
ing four-point mosaic tile pattern shown in Figure 12. 



energy (keV) energy (keV) 

Fig. 10.— The flux spectra generated using the APEC thermal 
plasma code (solid) for cluster 2 at 2 = 0.222 are shown along with 
the BAPEC fits to the Poisson spectra of Figure 9 (dashed). The 
spectra are generated with the cluster at a fiducial redshift of 0.1. 

The advantage of using simulations is that we can ex¬ 
amine the underlying cause of this observed velocity struc¬ 
ture. Figure 11 shows the density-weighted, LOS velocity 
map of the gas in cluster 2 in a region 2 r 2 oo on a side. 
A strong vertical gradient in projected velocity is appar¬ 
ent. The lower row of images in Figures 12 and 13 shows 
the z-axis (orthogonal view) surface brightness and tem¬ 


perature maps of cluster 2 at output redshifts z = 0.290, 
0.222 and 0.160. With this additional information, we can 
see that cluster 2 is extremely active. In this case it ap¬ 
pears that the redshifted quadrants (the lower half of the 
mosaic) represent the primary cluster itself, which has a 
peculiar velocity away from the observer. The blueshifted 
quadrants appear to contain a merging subcluster which is 
moving through the core in the negative y-direction. This 
subcluster can be traced in the z-axis projections starting 
in the upper left corner and moving toward the lower right 
over the course of the outputs presented in Figure 12. 



away from the observer. The map size is 2r2oo ^ side. 

3.2. Line Broadening 

To evaluate the significance of the observed line broad- 
enings, we normalize recovered from BAPEC by the 
expected thermal velocity dispersion {ath = i/feT/mpe)- 
The resulting velocity dispersion distribution is presented 
in Figure 14. The distribution ignores ~ 27% of the fits 
that return zero for the velocity broadening parameter (see 
§2.5). The distribution of non-zero values is fit reasonably 
well by 

> x) « 0.73 M -k j . (4) 

Out of the whole sample, 902 cluster projections, or 49%, 
yield a BAPEC velocity dispersion larger than the ther¬ 
mal broadening of the cluster. As in the gradient case, the 
probability drops as {pv!■ 

As we did for the velocity gradient study, we can also 
consider the impact of the epoch of observation and of the 
relative size of the Astro-E2 FOV on the results of Fig¬ 
ure 14. To determine if our results are dependent on the 
epoch of observation, we again construct three subsamples 
grouped by simulation output redshift and generate the cu¬ 
mulative frequency of normalized velocity broadenings for 
each subsample. The results are shown in Figure 15. The 
K-S test returns an 89.5% likelihood that the high redshift 
and middle redshift subsamples follow the same distribu¬ 
tion, a 29.3% likelihood that the middle and low redshift 
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Fig. 12.— Maps show the log of the unabsorbed soft (0.5-2 keV) X-ray surface brightness (arbitrary units) for cluster 2 at output redshifts 
0.290 (left), 0.222 (middle) and 0.160 (right), projected along the y (top) and 2 (bottom) axes of the simulation volume. The width of the 
image box is scaled to 2 r2oo each epoch. The mock Astro-E2 fields of view are overlaid on the frame that displays a velocity gradient 
An = 1.47cs. Note that although this map places the observer along the positive y-axis, the spectra were generated with positive y-axis 
velocity yielding redshift and negative y-axis velocity giving blueshift. 



Fig. 13.— Temperature maps (keV) for cluster 2 at output redshifts 0.290 (left), 0.222 (middle), and 0.160 (right) along the x (top) and y 
(bottom) axes. 





fraction of sample fraction of sample 
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subsamples follow the same distribution, and a 26.2% like¬ 
lihood that the high and low redshift subsamples follow the 
same distribution. 



0.01 0.10 1.00 10.00 


Fig. 14.— Fraction of sample with normalized velocity broaden¬ 
ing greater than the x-axis value. The dashed line represents the fit 
in equation (4). 



0.01 0.10 1.00 10.00 

Oy/Gth 


Fig. 15.— Cumulative likelihoods of the normalized velocity 
broadening for three subsamples grouped according to simulation 
output redshift. The short-dashed line is the subsample of clusters 
at ^ > 0.365, the long-dashed line is clusters at 0.365 > 2 : > 0.102, 
and the solid line is clusters at 2 : < 0.102. 

To evaluate the significance of the size of the Astro-E2 
FOV on the observed clusters, we create two subsamples 
grouped by the physical size of the cluster. The result¬ 
ing cumulative frequency distributions are shown in Fig¬ 
ure 16. For these distributions a K-S test is inappropriate 
because the result is dominated by the different sizes of 
the “zero-broadening” component of each subsample. Re¬ 
call from §2.5 that clusters with velocity broadenings less 
than about 50 km s“^ are expected to be unresolved by the 


XRS instrument. Clusters with smaller physical sizes are 
expected to have smaller physical broadenings. This fact 
implies that the small cluster subsample should exhibit a 
larger fraction of broadenings below the resolution of the 
XRS instrument than the large cluster subsample does. 
This effect is seen in the low-broadening limiting values 
of the distributions of Figure 16. This difference domi¬ 
nates the K-S test, resulting in a likelihood of 1.5 x 10“^ 
that the large cluster and small cluster subsamples follow 
the same distribution. If we perform the K-S test without 
including differences in the region Uy/ath < 2 we find a 
98% likelihood that the distributions are the same, imply¬ 
ing that the apparent divergence at large broadenings is 
statistically insignificant. 



0.01 0.10 1.00 10.00 


Fig. 16.— Cumulative likelihoods of the normalized velocity 
broadening for two equal subsamples grouped according to cluster 
size. The dashed line is the subsample of clusters with r2oo < 1-46 
Mpc, the solid line is clusters with r200 > 1-46 Mpc. 

Once again, we select one cluster projection to highlight 
as a case study. We select the cluster which exhibited 
the most significant line broadening: cluster 47 at output 
redshift 0.222 viewed along the x-axis. The data relevant 
to this cluster is compiled in Figures 17-20. Again it is 
apparent that the cluster is undergoing a strong merger 
at this epoch. In this case, however, the merging cluster 
is almost directly behind cluster 47 (rather than slightly 
off-center as in the case of cluster 2 output redshift 0.222 
shown above). This on-axis configuration means that all 
Astro-E2 spectra in the four-pointing mosaic contain emis¬ 
sion from the merging cluster, resulting in nearly identical 
(two-component) spectra. For this reason, each image in 
the multiple-pointing study gives approximately the same 
redshift when fit with a single-component plasma emission 
model, and differencing multiple images yields small veloc¬ 
ity gradients. In fact, while this projection gives a velocity 
dispersion more than ten times the thermal velocity broad¬ 
ening, its maximum velocity gradient in the four-pointing 
mosaic is 0.63cs. This illustrates the complementarity of 
line broadening and multiple-imaging studies for detecting 
dynamically active clusters. 
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6.00 6.10 6.20 6.30 

energy (keV) 


Fig. 17.— From top to bottom: the integrated spectrum, BAPEC 
fit, and actual flux spectrum for cluster 47, output at redshift 0.222, 
taken along the x-axis. The spectra are generated with the cluster 
at a fiducial redshift of 0.1. 


see file figure 18 .jpg 


Fig. 18.— Radial velocity map for cluster 47, output redshift 
0.222 projected along the x-axis in the same format as Figure 11. 

4. CONCLUSIONS 

We use a large statistical sample of mock Astro-E2 ob¬ 
servations to investigate the expected frequency of detec¬ 
tion of bulk motion of the ICM using the Fe XXV Ka line. 
We analyze differences in mean velocities among adjoining 


pointings as well as velocity broadenings in a single cen¬ 
tral pointing of clusters developed in a ACDM model and 
scaled to an observed redshift 2 : = 0.1. The mean velocities 
and broadenings are recovered from BAPEC fits to mock 
Astro-E2 spectra assuming 200 counts and 400 counts in 
the line region, respectively. 

From the central pointings, we find that roughly one- 
half of randomly selected clusters display velocity broad¬ 
ening in the Fe Ka line at least as large as the thermal 
width, and more than 15% will be broadened at a level 
greater than 2ath- In the four-point mosaic case, 6% of 
clusters show a maximum velocity difference that is tran¬ 
sonic {Av/cs > 0.5) while only 0.7% will be supersonic 
{Av/cs > 1 ). 

These fractions were obtained without any attempt to 
select for interesting clusters or to optimally align obser¬ 
vations on the basis of morphology. Further, our sample 
is biased toward cooler (smaller) clusters. We have at¬ 
tempted to remove the effects of this bias by normaliza¬ 
tion, but the physical cut represented by the ~ 75 km s“^ 
resolution of the XRS instrument cannot be normalized 
away. We can only assume that an observational survey, 
conducted to optimize results by selecting hot clusters with 
bimodal galaxy velocity distributions would obtain a sam¬ 
ple with an even larger fraction of dynamically interesting 
clusters. 

In extreme cases, however, morphological clues from X- 
ray observations can be misleading. We have illustrated 
one example in which the viewing angle is nearly aligned 
with the infall direction of a merging satellite. Such cases 
possess relatively simple X-ray morphologies but show 
large velocity effects in the Fe Ka line. 

This calibration exercise bears out the expectation that 
Astro-E2 will be an important instrument for detecting 
bulk motion of ICM gas. With reasonable counting statis¬ 
tics, the true emission-weighted velocity within the Astro- 
E2 FOV can be obtained to within ~ 50 km s“^, assuming 
the pre-launch gain and response estimates are accurate 
and that background is negligible. 
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FITS format spectra and maps similar to those pre¬ 
sented in this paper are available on the VGE website for 
selected cluster projections used in this study, including 
all those which exhibited splittings greater than 0.75cs or 
broadenings larger than 2.5ath- 
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Fig. 19.— Maps show the log of the unabsorbed soft (0.5-2 keV) X-ray surface brightness (arbitrary units) for cluster 47 at output redshifts 
0.290 (left), 0.222 (middle) and 0.160 (right), projected along the x (top) and y (bottom) axes of the simulation volume. The width of the 
image box is scaled to 2r2oo each epoch. The mock Astro-E2 fields of view are overlaid on the frame that displays the largest velocity 
broadening of the sample. Note that although this map places the observer along the positive x-axis, the spectra were generated with positive 
x-axis velocity yielding redshift and negative x-axis velocity giving blueshift. 



Fig. 20.— Temperature maps (keV) for cluster 47 at output redshifts 0.290 (left), 0.222 (middle) and 0.160 (right) along the x (top) and y 
(bottom) axes. 
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APPENDIX 


COMPUTING OPTICAL DEPTH 


The optical depth for line emission obeys (Spitzer 1978): 


J T{v)dv = fjkNj 


ire 


rrieC 


1 _ 


(Al) 


where j and k label atomic states, Vjk is the central frequency of the line radiation emitted in a transition from k j, fjk 
is the upward oscillator stregth for the transition, and Nj is the column density of atoms/ions in state j. The quantities 
bk and bj measure the particle densities Uk and rij relative to the values they would have in thermodynamic equilibrium 
(see Spitzer (1978) for details). At ICM densities, we expect bk = 0 for the Fe Ka line (Dere et al. 1997; Young et al. 
2003). The oscillator strength for the Ka line of Fe XXV is 0.798 (Verner et al. 1996). In cgs units, then: 


J T{v)dv « 0.0212iV, 


(A2) 


where Nj is the column density (in cm ^) of Fe XXV in the ground state. 

We next simplify things by defining an average optical depth for the line: 

2'/2fvjk-^= J T{v)dv (A3) 

where we define a Gaussian sigma parameterizing the width of the Fe Ka line. In the terminology of this paper, we use 
the definition: 

^jk = \j (A4) 


(note that, in keeping with the rest of this work, we have chosen a cr in velocity units so that the 
unitless). We have computed (jjk for each cluster as a part of our line broadening study (see §3.2) 

(again, in cgs units): 

_ _ 4.6 X IQ-^^cNj 

V^th + 

The final quantity remaining to compute is Nj. 

Assuming an iron abundance of 0.4 solar, we can write: 

iVFe « 2 X lO-'^iV^. (A6) 

Not all the iron will be in the FeXXV ionization state, however. At maximum, 78% of the Fe will be Helium-like. This 
fraction has a non-negligible temperature dependence over the range of temperatures expected in the ICM. The peak 
occurs at a temperature of 3 keV. For temperatures below 1 keV or above 8 keV the fraction has fallen below 33%. For 
this reason, we use our knowledge of the underlying simulation to approximate the optical depth for resonant scattering 
in our sample of clusters. We define a function /(T) using a linear interpolation on Table 4 of Arnaud & Raymond (1992). 
We then assign the ith gas parcel within the FOV a column density for resonant scattering given by: 

N^^exxv = 2 X 10-^fiT,)Nh. (A7) 


combination cr/c is 
We can now write 

(A5) 


We can then sum up the column densities in the FOV to find the average optical depth in the FOV: 


TFOV = 9.2 X 10 ^ ( 

ieFOV 


cfiT.)Nh 

+ <^v 


(AS) 
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